Guild and Niche Determination Enable Targeted Alteration of the Microbiome.

Microbiome science has greatly contributed to our understanding of microbial life and its essential roles for the environment and human health1-5. However, the nature of microbial interactions and how microbial communities respond to perturbations remains poorly understood, resulting in an often descriptive and correlation-based approach to microbiome research6-8. Achieving causal and predictive microbiome science would require direct functional measurements in complex communities to better understand the metabolic role of each member and its interactions with others. In this study we present a new approach that integrates transcription and translation measurements to predict competition and substrate preferences within microbial communities, consequently enabling the selective manipulation of the microbiome. By performing metatranscriptomic (metaRNA-Seq) and metatranslatomic (metaRibo-Seq) analysis in complex samples, we classified microbes into functional groups (i.e. guilds) and demonstrated that members of the same guild are competitors. Furthermore, we predicted preferred substrates based on importer proteins, which specifically benefited selected microbes in the community (i.e. their niche) and simultaneously impaired their competitors. We demonstrated the scalability of microbial guild and niche determination to natural samples and its ability to successfully manipulate microorganisms in complex microbiomes. Thus, the approach enhances the design of pre- and probiotic interventions to selectively alter members within microbial communities, advances our understanding of microbial interactions, and paves the way for establishing causality in microbiome science.

A Genome-Scale Atlas Reveals Complex Interplay of Transcription and Translation in an Archaeon.

The scale of post-transcriptional regulation and the implications of its interplay with other forms of regulation in environmental acclimation are underexplored for organisms of the domain Archaea. Here, we have investigated the scale of post-transcriptional regulation in the extremely halophilic archaeon Halobacterium salinarum NRC-1 by integrating the transcriptome-wide locations of transcript processing sites (TPSs) and SmAP1 binding, the genome-wide locations of antisense RNAs (asRNAs), and the consequences of RNase_2099C knockout on the differential expression of all genes. This integrated analysis has discovered that 54% of all protein-coding genes in the genome of this haloarchaeon are likely targeted by multiple mechanisms for putative post-transcriptional processing and regulation, with about 20% of genes likely being regulated by combinatorial schemes involving SmAP1, asRNAs, and RNase_2099C. Comparative analysis of mRNA levels (transcriptome sequencing [RNA-Seq]) and protein levels (sequential window acquisition of all theoretical fragment ion spectra mass spectrometry [SWATH-MS]) for 2,579 genes over four phases of batch culture growth in complex medium generated additional evidence for the conditional post-transcriptional regulation of 7% of all protein-coding genes. We demonstrate that post-transcriptional regulation may act to fine-tune specialized and rapid acclimation to stressful environments, e.g., as a switch to turn on gas vesicle biogenesis to promote vertical relocation under anoxic conditions and modulate the frequency of transposition by insertion sequence (IS) elements of the IS200/IS605, IS4, and ISH3 families. Findings from this study are provided as an atlas in a public Web resource (https://halodata.systemsbiology.net). IMPORTANCE While the transcriptional regulation landscape of archaea has been extensively investigated, we currently have limited knowledge about post-transcriptional regulation and its driving mechanisms in this domain of life. In this study, we collected and integrated omics data from multiple sources and technologies to infer post-transcriptionally regulated genes and the putative mechanisms modulating their expression at the protein level in Halobacterium salinarum NRC-1. The results suggest that post-transcriptional regulation may drive environmental acclimation by regulating hallmark biological processes. To foster discoveries by other research groups interested in the topic, we extended our integrated data to the public in the form of an interactive atlas (https://halodata.systemsbiology.net).

synDNA-a Synthetic DNA Spike-in Method for Absolute Quantification of Shotgun Metagenomic Sequencing.

Microbiome studies have the common goal of determining which microbial taxa are present, respond to specific conditions, or promote phenotypic changes in the host. Most of these studies rely on relative abundance measurements to drive conclusions. Inherent limitations of relative values are the inability to determine whether an individual taxon is more or less abundant and the magnitude of this change between the two samples. These limitations can be overcome by using absolute abundance quantifications, which can allow for a more complete understanding of community dynamics by measuring variations in total microbial loads. Obtaining absolute abundance measurements is still technically challenging. Here, we developed synthetic DNA (synDNA) spike-ins that enable precise and cost-effective absolute quantification of microbiome data by adding defined amounts of synDNAs to the samples. We designed 10 synDNAs with the following features: 2,000-bp length, variable GC content (26, 36, 46, 56, or 66% GC), and negligible identity to sequences found in the NCBI database. Dilution pools were generated by mixing the 10 synDNAs at different concentrations. Shotgun metagenomic sequencing showed that the pools of synDNAs with different percentages of GC efficiently reproduced the serial dilution, showing high correlation (r = 0.96; R2 ≥ 0.94) and significance (P < 0.01). Furthermore, we demonstrated that the synDNAs can be used as DNA spike-ins to generate linear models and predict with high accuracy the absolute number of bacterial cells in complex microbial communities. IMPORTANCE The synDNAs designed in this study enable accurate and reproducible measurements of absolute amount and fold changes of bacterial species in complex microbial communities. The method proposed here is versatile and promising as it can be applied to bacterial communities or genomic features like genes and operons, in addition to being easily adaptable by other research groups at a low cost. We also made the synDNAs' sequences and the plasmids available to encourage future application of the proposed method in the study of microbial communities.

Zebra: Static and Dynamic Genome Cover Thresholds with Overlapping References.

Assigning taxonomy remains a challenging topic in microbiome studies, due largely to ambiguity of reads which overlap multiple reference genomes. With the Web of Life (WoL) reference database hosting 10,575 reference genomes and growing, the percentage of ambiguous reads will only increase. The resulting artifacts create both the illusion of co-occurrence and a long tail end of extraneous reference hits that confound interpretation. We introduce genome cover, the fraction of reference genome overlapped by reads, to distinguish these artifacts. We show how to dynamically predict genome cover by read count and examine our model in Staphylococcus aureus monoculture. Our modeling cleanly separates both S. aureus and true contaminants from the false artifacts of reference overlap. We next introduce saturated genome cover, the true fraction of a reference genome overlapped by sample contents. Genome cover may not saturate for low abundance or low prevalence bacteria. We assuage this worry with examination of a large human fecal data set. By compositing the metric across like samples, genome cover saturates even for rare species. We note that it is a threshold on saturated genome cover, not genome cover itself, which indicates a spurious reference hit or distant relative. We present Zebra, a method to compute and threshold the genome cover metric across like samples, a recurrence to estimate genome cover and confirm saturation, and provide guidance for choosing cover thresholds in real world scenarios. Standalone genome cover and integration into Woltka are available: https://github.com/biocore/zebra_filter, https://github.com/qiyunzhu/woltka. IMPORTANCE Taxonomic assignment, assigning sequences to specific taxonomic units, is a crucial processing step in microbiome analyses. Issues in taxonomic assignment affect interpretation of what microbes are present in each sample and may be associated with specific environmental or clinical conditions. Assigning importance to a particular taxon relies strongly on independence of assigned counts. The false inclusion of thousands of correlated taxa makes interpretation ambiguous, leading to underconstrained results which cannot be reproduced. The importance sometimes attached to implausible artifacts such as anthrax or bubonic plague is especially problematic. We show that the Zebra filter retrieves only the nearest relatives of sample contents enabling more reproducible and biologically plausible interpretation of metagenomic data.

Acetate reprograms gut microbiota during alcohol consumption.

Liver damage due to chronic alcohol use is among the most prevalent liver diseases. Alcohol consumption frequency is a strong factor of microbiota variance. Here we use isotope labeled [1-13C] ethanol, metagenomics, and metatranscriptomics in ethanol-feeding and intragastric mouse models to investigate the metabolic impacts of alcohol consumption on the gut microbiota. First, we show that although stable isotope labeled [1-13C] ethanol contributes to fatty acid pools in the liver, plasma, and cecum contents of mice, there is no evidence of ethanol metabolism by gut microbiota ex vivo under anaerobic conditions. Next, we observe through metatranscriptomics that the gut microbiota responds to ethanol-feeding by activating acetate dissimilation, not by metabolizing ethanol directly. We demonstrate that blood acetate concentrations are elevated during ethanol consumption. Finally, by increasing systemic acetate levels with glyceryl triacetate supplementation, we do not observe any impact on liver disease, but do induce similar gut microbiota alterations as chronic ethanol-feeding in mice. Our results show that ethanol is not directly metabolized by the gut microbiota, and changes in the gut microbiota linked to ethanol are a side effect of elevated acetate levels. De-trending for these acetate effects may be critical for understanding gut microbiota changes that cause alcohol-related liver disease.

The Ubiquitous Human Skin Commensal Staphylococcus hominis Protects against Opportunistic Pathogens.

Staphylococcus hominis is frequently isolated from human skin, and we hypothesize that it may protect the cutaneous barrier from opportunistic pathogens. We determined that S. hominis makes six unique autoinducing peptide (AIP) signals that inhibit the major virulence factor accessory gene regulator (agr) quorum sensing system of Staphylococcus aureus. We solved and confirmed the structures of three novel AIP signals in conditioned medium by mass spectrometry and then validated synthetic AIP activity against all S. aureus agr classes. Synthetic AIPs also inhibited the conserved agr system in a related species, Staphylococcus epidermidis. We determined the distribution of S. hominis agr types on healthy human skin and found S. hominis agr-I and agr-II were highly represented across subjects. Further, synthetic AIP-II was protective in vivo against S. aureus-associated dermonecrotic or epicutaneous injury. Together, these findings demonstrate that a ubiquitous colonizer of human skin has a fundamentally protective role against opportunistic damage. IMPORTANCE Human skin is home to a variety of commensal bacteria, including many species of coagulase-negative staphylococci (CoNS). While it is well established that the microbiota as a whole maintains skin homeostasis and excludes pathogens (i.e., colonization resistance), relatively little is known about the unique contributions of individual CoNS species to these interactions. Staphylococcus hominis is the second most frequently isolated CoNS from healthy skin, and there is emerging evidence to suggest that it may play an important role in excluding pathogens, including Staphylococcus aureus, from colonizing or infecting the skin. Here, we identified that S. hominis makes 6 unique peptide inhibitors of the S. aureus global virulence factor regulation system (agr). Additionally, we found that one of these peptides can prevent topical or necrotic S. aureus skin injury in a mouse model. Our results demonstrate a specific and broadly protective role for this ubiquitous, yet underappreciated skin commensal.

Intestinal α1-2-Fucosylation Contributes to Obesity and Steatohepatitis in Mice.

BACKGROUND & AIMS: Fucosyltransferase 2 (Fut2)-mediated intestinal α1- 2-fucosylation is important for host-microbe interactions and has been associated with several diseases, but its role in obesity and hepatic steatohepatitis is not known. The aim of this study was to investigate the role of Fut2 in a Western-style diet-induced mouse model of obesity and steatohepatitis. METHODS: Wild-type (WT) and Fut2-deficient littermate mice were used and features of the metabolic syndrome and steatohepatitis were assessed after 20 weeks of Western diet feeding. RESULTS: Intestinal α1-2-fucosylation was suppressed in WT mice after Western diet feeding, and supplementation of α1-2-fucosylated glycans exacerbated obesity and steatohepatitis in these mice. Fut2-deficient mice were protected from Western diet-induced features of obesity and steatohepatitis despite an increased caloric intake. These mice have increased energy expenditure and thermogenesis, as evidenced by a higher core body temperature. Protection from obesity and steatohepatitis associated with Fut2 deficiency is transmissible to WT mice via microbiota exchange; phenotypic differences between Western diet-fed WT and Fut2-deficient mice were reduced with antibiotic treatment. Fut2 deficiency attenuated diet-induced bile acid accumulation by altered relative abundance of bacterial enzyme 7-α-hydroxysteroid dehydrogenases metabolizing bile acids and by increased fecal excretion of secondary bile acids. This also was associated with increased intestinal farnesoid X receptor/fibroblast growth factor 15 signaling, which inhibits hepatic synthesis of bile acids. Dietary supplementation of α1-2-fucosylated glycans abrogates the protective effects of Fut2 deficiency. CONCLUSIONS: α1-2-fucosylation is an important host-derived regulator of intestinal microbiota and plays an important role for the pathogenesis of obesity and steatohepatitis in mice.

The sum is greater than the parts: exploiting microbial communities to achieve complex functions.

Multi-species microbial communities are ubiquitous in nature. The widespread prevalence of these communities is due to highly elaborated interactions among their members thereby accomplishing metabolic functions that are unattainable by individual members alone. Harnessing these communal capabilities is an emerging field in biotechnology. The rational intervention of microbial communities for the purpose of improved function has been facilitated in part by developments in multi-omics approaches, synthetic biology, and computational methods. Recent studies have demonstrated the benefits of rational interventions to human and animal health as well as agricultural productivity. Emergent technologies, such as in situ modification of complex microbial community and community metabolic modeling, represent an avenue to engineer sustainable microbial communities. In this opinion, we review relevant computational and experimental approaches to study and engineer microbial communities and discuss their potential for biotechnological applications.

A comparison of DNA/RNA extraction protocols for high-throughput sequencing of microbial communities.

One goal of microbial ecology researchers is to capture the maximum amount of information from all organisms in a sample. The recent COVID-19 pandemic, caused by the RNA virus SARS-CoV-2, has highlighted a gap in traditional DNA-based protocols, including the high-throughput methods the authors previously established as field standards. To enable simultaneous SARS-CoV-2 and microbial community profiling, the authors compared the relative performance of two total nucleic acid extraction protocols with the authors' previously benchmarked protocol. The authors included a diverse panel of environmental and host-associated sample types, including body sites commonly swabbed for COVID-19 testing. Here the authors present results comparing the cost, processing time, DNA and RNA yield, microbial community composition, limit of detection and well-to-well contamination between these protocols.

A comparison of DNA/RNA extraction protocols for high-throughput sequencing of microbial communities.

One goal among microbial ecology researchers is to capture the maximum amount of information from all organisms in a sample. The recent COVID-19 pandemic, caused by the RNA virus SARS-CoV-2, has highlighted a gap in traditional DNA-based protocols, including the high-throughput methods we previously established as field standards. To enable simultaneous SARS-CoV-2 and microbial community profiling, we compare the relative performance of two total nucleic acid extraction protocols and our previously benchmarked protocol. We included a diverse panel of environmental and host-associated sample types, including body sites commonly swabbed for COVID-19 testing. Here we present results comparing the cost, processing time, DNA and RNA yield, microbial community composition, limit of detection, and well-to-well contamination, between these protocols. Accession numbers: Raw sequence data were deposited at the European Nucleotide Archive (accession#: ERP124610) and raw and processed data are available at Qiita (Study ID: 12201). All processing and analysis code is available on GitHub ( github.com/justinshaffer/Extraction_test_MagMAX ). Methods summary: To allow for downstream applications involving RNA-based organisms such as SARS-CoV-2, we compared the two extraction protocols designed to extract DNA and RNA against our previously established protocol for extracting only DNA for microbial community analyses. Across 10 diverse sample types, one of the two protocols was equivalent or better than our established DNA-based protocol. Our conclusion is based on per-sample comparisons of DNA and RNA yield, the number of quality sequences generated, microbial community alpha- and beta-diversity and taxonomic composition, the limit of detection, and extent of well-to-well contamination.

Bacterial modification of the host glycosaminoglycan heparan sulfate modulates SARS-CoV-2 infectivity.

The human microbiota has a close relationship with human disease and it remodels components of the glycocalyx including heparan sulfate (HS). Studies of the severe acute respiratory syndrome coronavirus (SARS-CoV-2) spike protein receptor binding domain suggest that infection requires binding to HS and angiotensin converting enzyme 2 (ACE2) in a codependent manner. Here, we show that commensal host bacterial communities can modify HS and thereby modulate SARS-CoV-2 spike protein binding and that these communities change with host age and sex. Common human-associated commensal bacteria whose genomes encode HS-modifying enzymes were identified. The prevalence of these bacteria and the expression of key microbial glycosidases in bronchoalveolar lavage fluid (BALF) was lower in adult COVID-19 patients than in healthy controls. The presence of HS-modifying bacteria decreased with age in two large survey datasets, FINRISK 2002 and American Gut, revealing one possible mechanism for the observed increase in COVID-19 susceptibility with age. In vitro , bacterial glycosidases from unpurified culture media supernatants fully blocked SARS-CoV-2 spike binding to human H1299 protein lung adenocarcinoma cells. HS-modifying bacteria in human microbial communities may regulate viral adhesion, and loss of these commensals could predispose individuals to infection. Understanding the impact of shifts in microbial community composition and bacterial lyases on SARS-CoV-2 infection may lead to new therapeutics and diagnosis of susceptibility.

Transcriptional profiling of lung macrophages identifies a predictive signature for inflammatory lung disease in preterm infants.

Lung macrophages mature after birth, placing newborn infants, particularly those born preterm, within a unique window of susceptibility to disease. We hypothesized that in preterm infants, lung macrophage immaturity contributes to the development of bronchopulmonary dysplasia (BPD), the most common serious complication of prematurity. By measuring changes in lung macrophage gene expression in preterm patients at risk of BPD, we show here that patients eventually developing BPD had higher inflammatory mediator expression even on the first day of life. Surprisingly, the ex vivo response to LPS was similar across all samples. Our analysis did however uncover macrophage signature genes whose expression increased in the first week of life specifically in patients resilient to disease. We propose that these changes describe the dynamics of human lung macrophage differentiation. Our study therefore provides new mechanistic insight into both neonatal lung disease and human developmental immunology.

Genomic and Transcriptomic Evidence Supports Methane Metabolism in Archaeoglobi.

Euryarchaeal lineages have been believed to have a methanogenic last common ancestor. However, members of euryarchaeal Archaeoglobi have long been considered nonmethanogenic and their evolutionary history remains elusive. Here, three high-quality metagenomic-assembled genomes (MAGs) retrieved from high-temperature oil reservoir and hot springs, together with three newly assembled Archaeoglobi MAGs from previously reported hot spring metagenomes, are demonstrated to represent a novel genus of Archaeoglobaceae, "Candidatus Methanomixophus." All "Ca Methanomixophus" MAGs encode an M methyltransferase (MTR) complex and a traditional type of methyl-coenzyme M reductase (MCR) complex, which is different from the divergent MCR complexes found in "Ca Polytropus marinifundus." In addition, "Ca Methanomixophus dualitatem" MAGs preserve the genomic capacity for dissimilatory sulfate reduction. Comparative phylogenetic analysis supports a laterally transferred origin for an MCR complex and vertical heritage of the MTR complex in this lineage. Metatranscriptomic analysis revealed concomitant in situ activity of hydrogen-dependent methylotrophic methanogenesis and heterotrophic fermentation within populations of "Ca Methanomixophus hydrogenotrophicum" in a high-temperature oil reservoir.IMPORTANCE Current understanding of the diversity, biology, and ecology of Archaea is very limited, especially considering how few of the known phyla have been cultured or genomically explored. The reconstruction of "Ca Methanomixophus" MAGs not only expands the known range of metabolic versatility of the members of Archaeoglobi but also suggests that the phylogenetic distribution of MCR and MTR complexes is even wider than previously anticipated.

Interplay of Staphylococcal and Host Proteases Promotes Skin Barrier Disruption in Netherton Syndrome.

Netherton syndrome (NS) is a monogenic skin disease resulting from loss of function of lymphoepithelial Kazal-type-related protease inhibitor (LEKTI-1). In this study we examine if bacteria residing on the skin are influenced by the loss of LEKTI-1 and if interaction between this human gene and resident bacteria contributes to skin disease. Shotgun sequencing of the skin microbiome demonstrates that lesional skin of NS subjects is dominated by Staphylococcus aureus (S. aureus) and Staphylococcus epidermidis (S. epidermidis). Isolates of either species from NS subjects are able to induce skin inflammation and barrier damage on mice. These microbes promote skin inflammation in the setting of LEKTI-1 deficiency due to excess proteolytic activity promoted by S. aureus phenol-soluble modulin α as well as increased bacterial proteases staphopain A and B from S. aureus or EcpA from S. epidermidis. These findings demonstrate the critical need for maintaining homeostasis of host and microbial proteases to prevent a human skin disease.

Author Correction: Environmental stimuli drive a transition from cooperation to competition in synthetic phototrophic communities.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Environmental stimuli drive a transition from cooperation to competition in synthetic phototrophic communities.

Phototrophic communities of photosynthetic algae or cyanobacteria and heterotrophic bacteria or fungi are pervasive throughout the environment1. How interactions between members contribute to the resilience and affect the fitness of phototrophic communities is not fully understood2,3. Here, we integrated metatranscriptomics, metabolomics and phenotyping with computational modelling to reveal condition-dependent secretion and cross-feeding of metabolites in a synthetic community. We discovered that interactions between members are highly dynamic and are driven by the availability of organic and inorganic nutrients. Environmental factors, such as ammonia concentration, influenced community stability by shifting members from collaborating to competing. Furthermore, overall fitness was dependent on genotype and streamlined genomes improved growth of the entire community. Our mechanistic framework provides insights into the physiology and metabolic response to environmental and genetic perturbation of these ubiquitous microbial associations.

Gut bacteria responding to dietary change encode sialidases that exhibit preference for red meat-associated carbohydrates.

Dietary habits have been associated with alterations of the human gut resident microorganisms contributing to obesity, diabetes and cancer1. In Western diets, red meat is a frequently eaten food2, but long-term consumption has been associated with increased risk of disease3,4. Red meat is enriched in N-glycolylneuraminic acid (Neu5Gc) that cannot be synthesized by humans5. However, consumption can cause Neu5Gc incorporation into cell surface glycans6, especially in carcinomas4,7. As a consequence, an inflammatory response is triggered when Neu5Gc-containing glycans encounter circulating anti-Neu5Gc antibodies8,9. Although bacteria can use free sialic acids as a nutrient source10-12, it is currently unknown if gut microorganisms contribute to releasing Neu5Gc from food. We found that a Neu5Gc-rich diet induces changes in the gut microbiota, with Bacteroidales and Clostridiales responding the most. Genome assembling of mouse and human shotgun metagenomic sequencing identified bacterial sialidases with previously unobserved substrate preference for Neu5Gc-containing glycans. X-ray crystallography revealed key amino acids potentially contributing to substrate preference. Additionally, we verified that mouse and human sialidases were able to release Neu5Gc from red meat. The release of Neu5Gc from red meat using bacterial sialidases could reduce the risk of inflammatory diseases associated with red meat consumption, including colorectal cancer4 and atherosclerosis13.

Establishing microbial composition measurement standards with reference frames.

Differential abundance analysis is controversial throughout microbiome research. Gold standard approaches require laborious measurements of total microbial load, or absolute number of microorganisms, to accurately determine taxonomic shifts. Therefore, most studies rely on relative abundance data. Here, we demonstrate common pitfalls in comparing relative abundance across samples and identify two solutions that reveal microbial changes without the need to estimate total microbial load. We define the notion of "reference frames", which provide deep intuition about the compositional nature of microbiome data. In an oral time series experiment, reference frames alleviate false positives and produce consistent results on both raw and cell-count normalized data. Furthermore, reference frames identify consistent, differentially abundant microbes previously undetected in two independent published datasets from subjects with atopic dermatitis. These methods allow reassessment of published relative abundance data to reveal reproducible microbial changes from standard sequencing output without the need for new assays.

The impact of skin care products on skin chemistry and microbiome dynamics.

BACKGROUND: Use of skin personal care products on a regular basis is nearly ubiquitous, but their effects on molecular and microbial diversity of the skin are unknown. We evaluated the impact of four beauty products (a facial lotion, a moisturizer, a foot powder, and a deodorant) on 11 volunteers over 9 weeks. RESULTS: Mass spectrometry and 16S rRNA inventories of the skin revealed decreases in chemical as well as in bacterial and archaeal diversity on halting deodorant use. Specific compounds from beauty products used before the study remain detectable with half-lives of 0.5-1.9 weeks. The deodorant and foot powder increased molecular, bacterial, and archaeal diversity, while arm and face lotions had little effect on bacterial and archaeal but increased chemical diversity. Personal care product effects last for weeks and produce highly individualized responses, including alterations in steroid and pheromone levels and in bacterial and archaeal ecosystem structure and dynamics. CONCLUSIONS: These findings may lead to next-generation precision beauty products and therapies for skin disorders.